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Heading date is a critical trait for adaptation of rice to different cultivation areas and cropping seasons. We 
evaluated the heading dates of 1,123 chromosome segments substitution lines (CSSLs) in the genetic 
background of an elite rice variety Huajingxian74 (HJX74). A CSSL with the substituted segments from 
ZihuilOO exhibited late heading under both natural long-day (NLD) and natural short-day (NSD) 
conditions, and the late heading phenotype was controlled by two novel epistatic loci on chromosome 8 and 
chromosome 3, respectively, termed LH8 and EH3. The function of EH3 was dependent on the LH8 
genotype through epistatic interaction between EH3 Zihui100 and LH8 Zihui100 alleles. Genetic and molecular 
characterization revealed LH8 encodes a CCAAT-box-binding transcription factor with Heading datel 
(Hdl)-binding activity and may delay flowering by repressing the expression of Early heading datel (Ehdl). 
Our work provides a solid foundation for further study on gene interaction in heading date and has 
application in breeding rice with greater adaptability. 



Rice is a short- day (SD) plant, and is becoming an increasingly important model monocot plant for molecu- 
lar biological study on agronomic important traits such as heading date. Heading date, flowering time in 
cereal crops, is a critical agronomical determinant for adaptation to specific cropping locations and growing 
seasons for current varieties of cultivated rice and is an important consideration for rice breeders. The molecular 
and developmental determinants of heading date have thus particularly important genetic targets for domest- 
ication and for breeding new varieties of rice. 

The timing of flower initiation is controlled by both environmental and internal signals 1,2 . By taking advantage 
of pioneering works in long-day (LD) plant Arabidopsis, many molecular components of rice flowering process 
initiated from photoperiod and light signals have been identified through natural variations, quantitative trait 
locus (QTL) analyses and by screening mutant populations (as summarized and reviewed in 3 ). Two independent 
flowering pathways, the conserved Hdl -dependent, and unique Ehdl -dependent pathways, control heading time 
in rice 4-7 . The two pathways function redundantly under SD conditions but antagonistically under LD condi- 
tions 8 . Hdl, a rice ortholog of Arabidopsis CONSTANS (CO) that is regulated by OsGI (a rice ortholog of 
Arabidopsis GIGANTEA), activates flowering under SD conditions and represses it strongly under LD conditions 
through regulating the expression of Heading date3a (Hd3a), a rice ortholog of Arabidopsis florigen 
FLOWERING LOCUS T (FT) 4 ' 9 . Ehdl encodes a B-type response regulator that is highly conserved in cultivated 
rice, but has no homolog in Arabidopsis 6 . Ehdl can promote flowering by activating the expression of Hd3a and 
RICE FLOWERING LOCUS Tl (RFT1), the closest paralog of Hd3a 6,7,10 . The expression of Ehdl is promoted by a 
number of positive regulators, including RID1 (also known as OsIDl and Ehd2) n ~ 13 , Early heading date3 (Ehd3) 14 , 
Early heading date4 (Ehd4) 15 , OsMADS51 16 and OsMADS50 (also known as DTH3) 17 ~ 19 . On the other hand, 
several factors repress the expression of Ehdl, thereby delaying flowering, such as Ghd7 20 , DTH8 (also known as 
Ghd8 and LHDl) 2123 , SE5 24 , OsCOL4 25 and OsLFLl 2627 . Therefore, Ehdl serves as a critical convergence point of 
regulation by multiple signaling pathways. The flowering signals are collected from Hdl and Ehdl, and are 
transducted to a pair of florigen genes, Hd3a and RFT1, promoting heading under SD and LD conditions, 
respectively 7,10 . After production in leaves, these florigens move to the shoot apical meristem and promote floral 
transition by activating the expression of OsMADS14 and OsMADS15 7 . Although these studies have revealed 
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much insight into the flowering control of rice, the underlying 
molecular mechanisms, especially the genetic interactions among 
these factors, are still not well understood. 

The advantages of CSSLs in precise and detailed phenotype evalu- 
ation, minor QTL detection and validation of gene-gene and gene- 
environment interactions facilitated both genetic studies and rice 
improvement 28 " 30 . Therefore, we had constructed a library of 1,123 
CSSLs in rice using HJX74, an elite indica variety from south China, 
as recipient and 24 accessions, including 14 indica and 10 japonica, 
collected worldwide as donors to understand the molecular mech- 
anism of agronomic traits including heading date 31,32 . Each CSSL 
contains a very small number of well -characterized chromosome 
segments from one of the 24 donor varieties. 

In this study, we found a CSSL with the substituted segments from 
ZihuilOO exhibited late heading under both NLD and NSD condi- 
tions, and the late heading phenotype was controlled by two novel 
epistatic loci on chromosome 8 and chromosome 3, respectively, 
termed LH8 and EH3. Map-based cloning of LH8 revealed it encodes 
a CCAAT-box-binding transcription factor and is allelic to DTH8. 
Molecular function of LH8 was further studied by expression analysis 
and yeast-two-hybrid experiment. 

Results 

Identification of a CSSL with late heading. The procedure for the 
development of the CSSLs was summarized in Fig. 1. HJX74 (as the 
female parent) was crossed with one of the 24 donors (ZihuilOO was 
taken as an example), and the F x plants were backcrossed with HJX74 
to develop the BQFi generation. Polymorphic SSR markers were 
used in the selection of the donor chromosomal segments. Using 
the same method, BC 6 Fi plants were obtained, and were self- 
crossed to produce BC 6 F 2 lines in which the majority of genomic 
regions were homozygous for HJX74 alleles. Totally, 1,123 such 
CSSLs in the HJX74 genetic background were developed 31 ' 32 . To 
further understand how flowering in rice was controlled, we 
evaluated the heading date of the 1,123 CSSLs under natural 
growth conditions. A CSSL, CSSL5, with substituted segments 
from ZihuilOO exhibited late heading compared with the recipient 
HJX74 under both NLD (114.4 ± 1.4 d for CSSL5, 100.1 ± 0.5 d for 



HJX74) and NSD conditions (84.5 ± 1.6 d for CSSL5, 76.7 ± 1.4 d for 
HJX74) in 2009 (Fig. 2a and 2b). 

To analyze the genetic basis for the late heading of CSSL5, we 
crossed CSSL5 with the recipient HJX74 (Fig. 1). 73 out of 326 F 2 
plants showed heading date later than that of CSSL5 (113.0 ± 1.3 d) 
under NLD conditions in 2010, and the whole population exhibited 
trimodal distribution of heading date fitting 4:9:3 segration ratio (x 2 
= 4.67 < x 2 o.o5,2 = 5.99) (Fig. 2c). These results indicate that the late 
heading phenotype was controlled by two genes with negative 
epistasis. 

Confirmation of the two genes by resequencing. To delimit the 
heading date genes on the substituted segments, CSSL5 and HJX74 
were high-throughput genotyped by whole-genome resequencing, 
and an ultrahigh- quality physical map was constructed (Fig. 3). 11 
substituted segments with a total length of 20.64 Mb derived from 
ZihuilOO were found to be distributed over 6 chromosomes (Fig. 3 
and Supplementary Table SI online). Molecular markers were 
subsequently designed on the substituted segments, and were used 
to analyze their linkage with heading date phenotype of the F 2 
population. Two QTLs were identified on chromosome 3 and 8, 
respectively, designated as EH3 and LH8, respectively. EH3 was 
mapped in the 750 kb genomic region between the Id32 and Id33 
markers and had the contributions to phenotypic variation by 14.0%, 
the allele from ZihuilOO could shorten the heading date by 2.32d and 
the dominant effect was —3.24. LH8 was mapped in the 460 kb 
genomic region between the Id83 and Id82 markers and had the 
contributions to phenotypic variation by 65.3%. The allele from 
ZihuilOO could delay the heading date by 6.62d and the dominant 
effect was 3.44 (Table 1). 

Epistatic interaction between EH3 and LH8. To examine the 
genetic interactions between EH3 and LH8, we scored heading 
dates in the F 2 population under NLD and NSD conditions in 
2011 (Fig. 4). The F 2 plants were classified into nine genotype 
classes based on the EH3 and LH8 alleles, and we compared the 
average heading date among the nine classes. Under NLD 
conditions, two-way ANOVA revealed the existence of digenic 
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Figure 1 | Flowchart of the development of CSSLs and QTL analysis. MAS: marker-assisted selection. 
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ZihuilOO allele at LH8 (i.e., increased days -to -heading (DTH)) was 
regardless of the allele status of EH3. In comparison, the effect of the 
ZihuilOO allele at EH3 (i.e., shortened DTH) was observed in two 
genotype classes, homozygous for the ZihuilOO allele at the LH8 
locus and heterozygous, but not in the class homozygous for the 
HJX74 allele at the LH8 locus (Fig. 4a). Under NSD conditions, the 
ZihuilOO allele at LH8 was also observed to delay heading regardless 
of the allele status of EH3. However, the two-way ANOVA suggested 
the digenic interaction between EH3 and LH8 was not significant 
under NSD conditions (P = 0.883) (Fig. 4b). 

Map-based cloning of LH8. To fine map the LH8 gene, we selected 
F 2 plants with homozygous EH3 allele from HJX74 and heterozygous 
LH8 allele from ZihuilOO to generate F 3 population (Fig. 1). F 3 
population showed bimodal distribution for DTH fitting 3 : 1 
segregation ratio under both NLD and NSD conditions in 2012 (% 2 
= 0.04 < x 2 o.()5,i = 3.84 and f = 0.03 < x 2 0 .o 5 ,i = 3.84 for NLD and 
NSD, respectively), indicating the ZihuilOO allele at LH8 increased 
DTH in a dominant manner (Fig. 5a and 5b). 7 new markers were 
developed in the marker interval of Id83-Id82, and were used to 
analyze a total of 2,159 F 3 plants. 6 recombinants were identified 
between markers RM22475 and RM25 (Fig. 5c). The self- 
pollinated progeny (F 4 lines) of those 6 plants were used to 
determine the genotypes of LH8 (Fig. 1). The recombinants r2 and 
r5 restricted LH8 to the 28.3 kb genomic region between the Id87 
and Id811 markers, including 3 putative genes, a transferase family 
protein (LOC_Os08g07730), a SERK-family receptor-like protein 
kinase (LOC_Os08g07760) and a CCAAT-box-binding transcrip- 
tion factor (LOC_Os08g07740) which had been reported to 
function in heading date control in rice (Fig. 5c). The sequencing 
of LH8 genomic regions from CSSL5 and HJX74 revealed two GGC 
insertion and a 1,116-bp deletion in LOC_Os08g07740 of HJX74 
(Fig. 6a), resulting in a 2-glycine insertion in the middle and a big 
alternation in the C-terminal region of the translated protein, 
respectively (Fig. 6b). These results suggest that LOC_Os08g07740 
corresponds to LH8 controlling the late heading in CSSL5. 
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Figure 2 | The heading dates of HJX74, CSSL5 and their F 2 population. 

(a) . Phenotypes of HJX74 (left) and CSSL5 (right). Yellow bar = 15 cm. 

(b) . Days-to-heading in HJX74 and CSSL5 in different photoperiodic 
conditions. NLD, natural long-day conditions; NSD, natural short-day 
conditions. Each column represents mean ± s.d. (n = 30). (c). Frenquency 
distributions of days-to-heading in an F 2 population from HJX74 x CSSL5 
under NLD conditions. Black arrows indicate the mean DTH of HJX74 and 
CSSL5 (n = 30), respectively. 

interaction between EH3 and LH8 (P = 5.96e-06). Orthogonal 
contrast test showed that A A (Additive by additive), but not AD 
(Additive by dominance), DA (Dominance by additive) or DD 
(Dominance by dominance) type of interaction was significant in 
the interaction of EH3 and LH8 (P < 0.0001). The effect of the 



Characterization of LH8 in rice. To examine whether the sequence 
mutations can alter the expression of LH8, we examined the 
expression levels of LH8 in CSSL5, HJX74 and NIL-LH8 (a line 
selected from LH8 F 2 segragating population with homozygous 
EH3 allele from HJX74 and homozygous LH8 allele from 
ZihuilOO) under NSD conditions in 2013. No significant difference 
in the LH8 expression levels could be detected among the three 
materials, suggesting that late heading did not result from the LH8 
transcription. The expression levels of other heading date genes were 
also analyzed. For Hdl and OsMADS50, no difference could be 
detected among the three materials under NSD conditions. The 
expression levels of Ehdl, Hd3a and RTF1 were less in NIL-LH8 
than those in HJX74 under NSD conditions, indicating that LH8 
may delay heading through suppressing the expressions of these 
floral transition activators. In contrast, the expression levels of 
Ehdl and RTF1 were higher in CSSL5 than those in HJX74. The 
different expression levels of Ehdl and RFT1 between CSSL5 and 
NIL-LH8 may result from the function of different EH3 alleles 
(Fig. 7a). 

To study the involvement of the 2-glycine indel and the C-ter- 
minal region of the LH8 protein in its function, we tested the inter- 
action between LH8 and Hdl by yeast-two-hybrid method. Yeast 
cells that coexpressed the LH8 Zihuil0 ° bait and Hdl prey fusion pro- 
teins were able to grow on media lacking histidine and adenine, 
indicating that these proteins can bind each other in yeast cells to 
promote the expression of the HIS3 and ADE2 reporter genes. 
However, the similar interaction could not be detected between 
LH8 HJX74 and Hdl, indicating the 2-glycine indel and the C-terminal 
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Figure 3 | Physical map of CSSL5 based on the whole-genome resequencing data. Red areas indicated the substituted chromosome segments from 
ZihuilOO, while the blue areas indicated chromosome regions of HJX74. The tiny yellow area on chromosome 9 indicated a heterozygous segment 
between HJX74 and ZihuilOO. 



Table 1 Heading date QTLs detected 
HJX74 x CSSL5 under NLD conditions 


in an F2 


population from 


Marker 

QTL Chromosome interval LOD a 


PVE b 


a c d d 


EH3 3 Id32-ld33 4.13 
LH8 8 Id83-ld82 53.12 


14.0% 
65.3% 


-2.32 -3.24 
6.62 3.44 


a Log-likelihood value. 

Percentage phenotypic variance explained by the QTL. 
c Additive effect of the Zihui 1 00 allele on days to heading. 
d Dominant effect of the Zihui 1 00 allele. 



region of LH8 might be important for its interaction with Hdl 
(Fig. 7b). 

Discussion 

Until now, a total of 734 rice heading date QTLs have been reported 
(http://www.gramene.org/qtl), and some of them have been molecu- 
larly characterized 3 . Besides understanding the biological function of 
single gene/QTL, clarification of genetic interactions among these 
genes/QTLs is also important. Actually, several examples of epistatic 
interactions among heading date genes/QTLs have been 
reported 33 " 37 . By eliminating genetic background noises, CSSLs offer 
a practical solution to bridge the huge gap of knowledge between the 
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Figure 4 | Epistatic interaction between LH8 and EH3. Differences in days-to-heading for nine genotype classes for combinations of LH8 and EH3 
under NLD (a) and NSD (b) conditions in the F 2 populations. Genotypes were determined using the closely linked markers Id83-Id82 (LH8) and Id32- 
Id33 (EH3). Z, Hetero and H indicate homozygous for ZihuilOO allele, heterozygous and homozygous for the HJX74 allele, respectively. Paa, 
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Figure 5 | Map-based cloning of the LH8 gene in rice, (a and b). Frenquency distributions of days-to-heading in an F 3 population from HJX74 x CSSL5 
under NLD (a) and NSD (b) conditions. Black arrows indicate the mean DTH of HJX74 and NIL-LH8 (n = 30), respectively, (c). High-resolution 
mapping of LH8. The genetic map of LH8 was based on recombinant events among 2,159 F 3 plants from HJX74 x CSSL5. The number of recombinants 
between adjacent markers is shown above the bar. White, shadow and black regions shown in the recombinants indicate homozygous regions for HJX74 
allele, heterozygous and homozygous regions for the ZihuilOO allele, respectively. 



genotype and phenotype, and therefore, are especially useful in 
studying complicated traits including heading date. In this study, a 
novel epistatic interaction between LH8 and EH3 involved in flower- 
ing control in rice was detected. Under NLD conditions and HJX74 
genetic background, in plants homozygous for the LH8 H]X74 allele, 
difference in DTH between EH3 ZihuiW0 and EH3 H]X74 alleles was small 
(about 1 day), as shown in HJX74 and NIL-EH3 (a line selected from 
LH8 F 2 segragating population with homozygous LH8 allele from 
HJX74 and homozygous EH3 allele from ZihuilOO) (Fig. 8a). A large 
effect of EH3 Zihui100 allele on heading promotion (about 7 days) 
occurred with the LH8 Zihuil0 ° allele, indicating EH3 Zihuil0 ° allele could 
function to shorten heading date only when LH8 Zihnil0 ° allele was 
concurrently present, and therefore, the function of EH3 was 
dependent on the LH8 genotype through epistatic interaction 
between EH3 Zihuil0 ° and LH8 Zihail0 ° alleles (Fig. 8a). This result sug- 
gests that EH3 might function as a modifier of LH8. Genetic inter- 
action between LH8 and EH3 suggested that QTL pyramiding would 
be an effective method for the development of varieties with desirable 
heading dates under different growth conditions. 

To facilitate the identification of the location of LH8 and EH3, the 
high-throughput genotyping technology was exploited. The recipi- 



ent HJX74 and CSSL5 were resequenced at the sequencing depth of 
66- and 1-fold, respectively. Besides a large substituted segment from 
the short arm end to RM282 on chromosome 3 which was detected 
through traditional marker- assisted selection during CSSLs con- 
struction 31 , 10 additional substituted segments were identified dis- 
tributing over 6 chromosomes, including some small introgressed 
segments that might be difficult to be detected by SSR markers (see 
Supplementary Table SI online). Indeed, LH8 was found on untar- 
geted chromosomal region on chromosome 8. These results showed 
that the high-throughput genotyping method by resequencing has a 
higher accuracy than marker- assisted selection. Owning to the deep 
sequencing coverage of the recipient HJX74, resequencing the CSSL5 
with genome coverage as low as 1-fold was enough to rapidly validate 
and delimit the QTL, suggesting the efficiency of this mapping-by- 
sequencing method in discovering the genes responsible for quant- 
itative trait loci. Our results presented a good example of combining 
the advanced material with uniform genetic background and high- 
throughput resequencing technology in economically accelerating 
gene identification. 

The sequencing of LH8 genomic regions from CSSL5 and HJX74 
revealed two GGC insertion and a 1,116-bp deletion in 
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HJX74 1 MKSRKSYGHLLSPVGSPPSDNESGAAAAAAAAGGGGCGSSAGYVVYGGGGGGDSPAKEQD 60 
CSSL5 1 MKSRKSYGHLLSPVGSPPSDNESGAAAAAAAAGGGGCGSSAGYVVYGGGGGGDSPAKEQD 60 

HJX74 61 RFLPIANVSRIMKRSLPANAKISKEAKETVQECVSEFISFVTGEASDKCQREKRKTINGD 120 
CSSL5 61 RFLPIANVSRIMKRSLPANAKISKEAKETVQECVSEFISFVTGEASDKCQREKRKTINGD 120 

HJX74 121 DLLWAMTTLGFEAYVGPLKSYLNRYREAEGEKAAVLGGAGGAAAARHGEGGCCGGGGGGG 180 
CSSL5 121 DLLWAMTTLGFEAYVGPLKSYLNRYREAEGEKAAVLGGAGGAAAARHGEGGCCGGGGG — 178 

HJX74 181 ADGVVIDGHSPLAGGLSHSHHGHQQQDGGGDVGLMMGGGAAGVGYNAGAGSTTTAFYAPA 240 
CSSL5 179 ADGVVIDGHSPLAGGLSHSHHGHQQQDGGGDVGLMMGGGAAGVGYNAGAGSTTTAFYAPA 238 

HJX74 241 ATAASGNKAYCGGDGSRVMEFEGIIVALYT AWRMLLYCPFSLPLLWLHA 289 

CSSL5 239 ATAASGNKAYCGGDGSRVMEFEGIGGEEESGGGGGGGERGFAGHLHGVQWFRLKRSTN 296 

Figure 6 | Sequence comparision of LH8 from HJX74 and CSSL5. (a). Structure of LH8 in HJX74 and CSSL5. Black and gray regions represent the ORF 
and UTR regions of LH8, respectively. The dotted line represents the 1,116 bp deletion in HJX74. (b). Protein alignment of LH8 from the predicted 
protein of HJX74 and CSSL5 using the EMBL software ClustalW2 multiple sequence alignment tool. 



LOC_Os08g07740 of HJX74, resulting in a 2-glycine insertion in the 
middle and a big alternation in the C-terminal region of the trans- 
lated protein, respectively (Fig. 6). Such sequence variations in the 
coding region of LOC_Os08g07740 between CSSL5 and HJX74 sug- 
gested that LOC_Os08g07740 might correspond to LH8 controlling 
the late heading in CSSL5 and NIL-LH8. LOC_Os08g07740 encodes 
a CCAAT-box-binding transcription factor which had been reported 
to strongly inhibit floral transition with names of DTH8> Ghd8 and 
LHDl 2123 . Similar to DTH8 and Ghd8, LH8 increased plant height 
and grain number per panicle under NLD conditions (Fig. 2 and 
Supplementary Table S2 online), indicating LH8 is allelic to DTH8, 
Ghd8 and LHDl. Indeed, the two GGC insertion and one 1,116-bp 
deletion in the coding region of LOC_Os08g07740 of HJX74 had also 
been found in the ZS97 allele at Ghd8 22 . However, unlike DTH8 and 
LHDl delaying flowering only under NLD/LD conditions 21,23 and 
Ghd8 functioning oppositely under SD and LD conditions 22 , LH8 
delays the heading date under both NSD and NLD conditions. 
These differences in function may come from different gene-back- 
ground interactions or different growth conditions. Therefore, genes 
controlling important agronomic traits are necessary to be further 
studied in different genetic backgrounds and under different natural 
field conditions to give a full understanding of its effects and to offer 
practical instructions in rice breeding. 

Expression analysis revealed the transcripts of several flowering 
promoting factors including Ehdl, Hd3a and RFT1 were fewer in 
NIL-LH8 than those in HJX74 under NSD conditions, suggesting 
LH8 might delay flowering by down-regulating the expression of 
Ehdl, Hd3a and RFT1 (Fig. 8b). The expressions levels of Hdl and 
OsMADS50 showed no significant differences among HJX74, CSSL5 
and NIL-LH8 under NSD conditions, suggesting they might function 
upstream of LH8 or parallel to LH8. Furthermore, no significant 
difference of LH8 expression were found among HJX74, CSSL5 
and NIL-LH8 under NSD conditions, indicating the reason for the 
change of heading date was not caused by LH8 expression, but by 
nucleotide changes of the coding sequence. The 6-bp and 1,116-bp 



indels in the coding region of LH8 resulted in a big sequence alterna- 
tion of the translated protein, especially in the C-terminal region. 

The protein encoded by LH8 is a putative HAP3 subunit of the 
HAP complex, which binds to the CCA AT box, a cis- acting element 
present in approximately 25% of eukaryotic gene promoters 38 ' 39 . 
HAP proteins always form a HAP2/HAP3/HAP5 trimeric complex 
with DNA binding activity in mammals 40,41 . In Arabidopsis, CO 
could replace AtHAP2 in the HAP complex to form a trimetric 
complex, CO/AtHAP3/AtHAP5, as the CCT domain of CO and 
the DNA-binding domain of AtHAP2 are similar to each other 42 . 
In LD conditions, overexpression of AtHAP2 oxAtHAP3 could delay 
flowering by impairing formation of the CO/AtHAP3/AtHAP5 
complex, leading to decreased expression of FT* 2 . As a homologue 
of CO, Hdl may be able to replace OsHAP2 in the HAP complex in 
rice. Our yeast-two-hybrid results revealed that Hdl HIX74 could 
indeed interact with LH8 Zihuil0 °, but not with LH8 HJX74 , indicating 
the sequence difference, especially in the C-terminal region, between 
LH8 zihuiioo and LH8 hjx74 might be i mportant f or i ts interaction with 

Hdl HJX74 , and therefore for the regulation of rice heading (Fig. 8b). 
The domains responsible for the interaction between LH8 HIX74 and 
Hdl HJX74 were needed to be further investigated. However, the inter- 
action between LH8 Zihuil0 ° and Hdl Nip (Hdl from Nipponbare) could 
not be detected (data not shown), in consistent with the previous 
study that effects of Ghd8 depend on the genetic background 22 . 
Consequently, whether the HAP complex could be formed or not 
was determined by the different subunits presented in different gen- 
etic backgrounds, possibly resulting in varied functions in flowering 
control. 

Until now, 4 genes (Ehd4, OsMADS50, OsDofl2, OsPhyB) on the 
short arm of chromosome 3 had been identified to be involved in 
flowering control in r i ce 15 > 17 - 19 > 43 > 44 . In our study, EH3 was mapped in 
the 750 kb genomic region between the Id32 and Id33 markers with 
123 putative ORFs inside, including OsMADS50. OsMADS50, an 
ortholog of Arabidopsis SUPPRESSOR OF OVEREXPRSSION OF 
CONSTANS1 (SOC1), functions upstream of Ehdl to indirectly 
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Figure 7 | Characterization of the LH8 gene in rice. (a). Expression of 
LH8, Hdl, OsMADS50, Ehdl, Hd3a and RFTl in HJX74 (white pillars), 
CSSL5 (gray pillars) and NIL-LH8 (black pillars) at dawn (for LH8, 
OsMADS50, Ehdl, Hd3a and RFTl) and dusk (for Hdl) under NSD 
conditions. The analysis was performed in two independent experiments. 
The mean values of the relative expression levels of genes in CSSL5 and 
NIL-LH8 compared with HJX74 are shown by cloumns, and standard 
deviations are indicated by the error bars. (b). Interaction of LH8 and Hdl 
by yeast- two -hybrid assay. The yeast cells were grown on SD/-Leu/-Trp/- 
His/-Ade medium. The interaction between SE and HYL1 was served as 
positive control, while the combination of empty prey pGADT7, and that 
of the empty bait vector pGBKT7 with/ without LH8 from ZihuilOO were 
served as negative controls. 

induce .RFTl transcription 17 " 19 and promote flowering. However, no 
epistatic interaction had been found between LH8 and OsMADS50. 
Furthermore, no significant difference of the expression of 
OsMADS50 was found between HJX74 and CSSL5. Therefore, 



whether EH3 is allelic to OsMADS50 remains to be investigated. 
Further fine mapping and molecular identification of EH3 will be 
necessary to provide molecular evidence for the epistatic interaction 
between LH8 and EH3 in regulating heading date in rice. 

Methods 

Plant materials and growth conditions. CSSL5 is one of a set of 1,123 chromosome 
segment substitution lines, which were developed from backcross progenies (BC 6 F 2 ) 
derived from a cross between an elite indica variety from south china, HJX74, as the 
recurrent parent and another indica variety ZihuilOO as the donor parent. 

Plants were grown under natural short-day (from July to November) and natural 
long-day (from March to July) conditions in a paddy field in Guangzhou (23°07'N, 
113°15'E), China. Heading date was defined as the time when the first panicle 
emerged. 

High-throughput genotyping using whole-genome resequencing. DNA was 

extracted from 100 mg fresh rice leaves using the DNeasy Plant Mini Kit (QIAGEN 
Sciences). DNA was quantified using NanoDrop 1000 spectrophotometer (Thermo 
Scientific), and about 5 ug of DNA was used for preparation of libraries for Illumina 
sequencing according to the protocol for the Paired-End DNA Sample Prep kit 
(Illumina). The libraries were used for cluster generation on a flow cell and sequenced 
for 76 cycles on an Illumina Genome Analyzer IIx. Base calling and filtering of low- 
quality bases were done using sequence control software real-time analysis, Base 
calling (BCL) converter and the GERALD module (Illumina). 

DNA extraction and PCR amplification. Fresh leaves were collected at the seedling 
stage and then ground in liquid nitrogen. Microquantities of DNA were extracted 
from fresh leaves of each individual using a previously reported method 45 . 
Amplification was carried out on the program for the initial denaturing step with 
94°C for 3 min, followed by 35 cycles for 30 s at 94°C, 30 s at 55°C, 30 s at 72°C, with 
a final extension at 72°C for 5 min. PCR products were separated on 6% non- 
denaturing polyacrylamide gel and detected using the silver staining method 46 . 

Linkage map and QTL analysis. We constructed the linkage map using Mapmaker/ 
Exp3.0 47 . The Kosambi function was used to calculate genetic distances. QTL for 
heading date were estimated by Mapmaker/QTL l.l 48 . Putative QTL were identified 
in regions exceeding 3.0 LOD (log-likelihood value). 

Evaluation of genetic interactions. To analyze genetic interactions between LH8 and 
EH3, we classified a total of 359 and 225 F 2 plants grown under NLD and NSD 
conditions in 2011, respectively, into nine classes according to the genotypes of SSR 
markers tightly linked to the QTLs (Id83-Id82 interval for LH8, and Id32-Id33 
interval for EH3). We compared the mean DTH between the nine genotypes using 
two-way analysis of variance (ANOVA). Testing for each genetic interaction types 
(AA, AD, DA and DD) was done by orthogonal contrast test in STAR software 
(http://bbi.irri.org/) . 

Mapping for LH8. F 2 plants with homozygous EH3 allele from HJX74 and 
heterozygous LH8 allele from ZihuilOO were selected to generate F 3 population. A 
total of 559 and 1,600 F 3 plants were grown to fine map LH8 gene under NLD and 
NSD conditions in 2012, respectively. Each plant was genotyped using molecular 
markers and phenotyped by visual examination of the first panicle emergence. The F 4 
recombinant lines derived from the F 3 population were used for progeny testing. 
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RNA Extraction and QRT-PCR. Total RNA from 30d-old leaves were isolated using 
TRIZOL reagent (Invitrogen) following the manufacturer's instruction. First-strand 
cDNA was reverse transcribed from DNasel-treated RNA with oligo-dT as the primer 
using ReverTra Ace kit (Toyobo). Gene expression was measured by QRT-PCR using 
the ABI 7500 system (Life technologies). The QRT-PCR was carried out in a total 
volume of 20 ul containing IX SYBR Green Master Mix (Life technologies). We 
normalized the expression levels by using UBQ5 gene as internal control. Each set of 
experiments was repeated three times, and the relative standard curve quantification 
method was used to evaluate quantitative variation. The QRT-PCR procedure was 
conducted at 94°C for 3 min, followed by 40 cycles at 94°C for 30 s, 58°C for 30 s, and 
72°C for 30 s. The QRT-PCR primers were listed in Supplemental Table S3 online. 

Yeast two-hybrid assay. The cDNA fragments encoding the entire putative LH8 and 
Hdl proteins were PCR-amplified, verified by sequencing and cloned into Ndel- 
EcoRI restriction sites of the MATCHMAKER two -hybrid prey vector pGADT7 or 
bait vector pGBKT7 (Clontech), respectively. The bait and prey constructs were co- 
transformed into the yeast strain AH 109, and the transformed cells were plated on 
synthetic defined medium containing all essential amino acids except leucine and 
tryptophan (SD/-Leu/-Trp) medium and incubated at 30°C for 3 d and then diluted 
and applied onto SD/-Leu/-Trp/-His/-Ade medium and SD/-Leu/-Trp medium (as 
loading control) and cultured at 30°C for 2-3 days. The interaction between HYL1 
and SE was used as a positive control 49 . The primers used in vector construction were 
listed in Supplemental Table S3 online. 
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